    clear;clc;close all
    addpath('./util');

%% 载入数据
    LS = loadtiff('.\data\20200623_1_LS.tif');

%% 可视化原始体数据
    Visualization3D(LS, [0.325 0.325 1], linspace(0,0.1,65535))

%% beads中心点提取
    centroids_subpixel = detect_beads_LS_2(LS);

%% 投影视图
    figure, imshow(imadjust(max(volData, [], 3))),title('LS beads xy');
    for i = 1:size(centroids_subpixel,1)
        hold on, plot(centroids_subpixel(i,1), centroids_subpixel(i,2), 'o', 'MarkerSize', 5);
    end

    figure, imshow(imadjust(squeeze(max(volData, [], 2)))), title('LS beads yz');
    for i = 1:size(centroids_subpixel,1)
        hold on, plot(centroids_subpixel(i,3), centroids_subpixel(i,2), 'o', 'MarkerSize', 3);
    end

    %% 3D显示体数据及中心坐标叠加
    Visualization3D(LS, [0.325 0.325 1], linspace(0,0.1,65535),1,centroids_subpixel,1)